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THEME 


The  Fluid  Dynamics  Panel  of  AGARD  arranged  a  Symposium  o..  “Applications  of  Computational  Fluid  Dynamics  in 
Aeronautics”,  on  7— 10  April  1 986  in  Aix-en-Provence,  France.  The  purpose  of  the  Symposium  was  to  provide  an 
assessment  of  the  status  of  CFD  in  aerodynamic  design  and  analysis,  with  an  emphasis  on  emerging  applications  of  advanced 
computational  techniques  to  complex  configurations.  Sessions  were  devoted  specifically  to  grid  generation,  methods  for 
inviscid  flows,  calculations  of  viscous- in  viscid  interactions,  and  methods  for  solving  the  Navier-Stokes  equations.  The  31 
papers  presented  at  the  meeting  are  published  in  AGARD  Conference  Proceedings  CP-412  and  are  listed  in  the  Appendix 
to  this  report.  A  brief  synopsis  of  each  paper  and  some  general  conclusions  and  recommendations  are  given  in  this 
evaluation  report. 


La  Commission  de  Dynamique  des  Fluides  de  l’  AG  ARD  a  organise  a  Aix  en  Provence  du  7  au  1 0  avril  1 986  un 
Symposium  sur  “Les  applications  des  calculs  de  dynamique  des  fluides  en  aeronautique".  Le  but  de  ce  symposium  etait  de 
faire  le  point  sur  les  possibility  des  methodes  de  calcul  pour  la  conception  et  la  prevision  en  aerodynamique  en  insistant 
particulierement  sur  les  applications  nouvelles  des  techniques  avancees  de  calcul  pour  les  configurations  complexes.  Des 
sessions  particulieres  ont  ete  consacree  a  la  construction  des  maillages,  aux  methodes  pour  les  ecoulements  non  visqueux, 
aux  calculs  des  interactions  fluide  parfait-fluide  visqueux  et  aux  methodes  pour  resoudre  les  equations  de  Navier-Stokes.  Les 
3 1  communications  presentees  au  congres  et  dont  la  liste  figure  en  annexe  au  rapport  ont  ete  publiees  par  TAG  ARD  en  un 
volume  qui  constitue  les  Actes  de  ce  symposium  et  porte  la  reference  CP-412.  Ce  rapport  devaluation  contient  une  breve 
analyse  de  chacune  des  communications  ainsi  que  les  conclusions  generales  et  les  recommandations  que  1’ensemble  de  ces 
contributions  a  inspirees. 
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TECHNICAL  EVALUATION  REPORT 

AGARO  Fluid  Dynamics  Panel  Symposium  on 
APPLICATIONS  of  COMPUTATIONAL  FLUID  DYNAMICS  In  AERONAUTICS 

W. J.  McCroskey 

U.S.  Army  Aerof llghtdynamlcs  Directorate 
NASA  Ames  Research  Center 
Moffett  Field*  California 


ABSTRACT 

The  Fluid  Dynamics  Panel  of  AGARO  arranged  a  Symposium  on  "Applications  of  Computational  Fluid  Dynam¬ 
ics  In  Aeronautics,"  on  7-10  April  1986  In  Aix-en-Provence,  France.  The  purpose  of  the  Symposium  was  to 
provide  an  assessment  of  the  status  of  CFD  In  aerodynamic  design  and  analysis,  with  an  emphasis  on  emerg¬ 
ing  applications  of  advanced  computational  techniques  to  complex  configurations.  Sessions  were  devoted 
specifically  to  grid  generation,  methods  for  Invlscld  flows,  calculations  of  vlscous-lnvlscld  Interac¬ 
tions,  and  methods  for  solving  the  Navler-Stokes  equations.  The  31  papers  presented  at  the  meeting  are 
published  In  AGARD  Conference  Proceedings  CP-412  and  are  listed  In  the  Appendix  of  this  report.  A  brief 
synopsis  of  each  paper  and  some  general  conclusions  and  recomendatlons  are  given  In  this  evaluation 
report. 

1.  INTRODUCTION 

The  58th  Meeting  of  the  AGARD  Fluid  Dynamics  Panel  (FDP)  was  held  from  the  7th  to  the  11th  of  April, 
1986,  In  Aix-en-Provence,  France.  It  Included  a  major  and  timely  symposium,  organized  by  the  recently 
established  FDP  standing  committee  on  Computational  Fluid  Dynamics,  with  the  following  theme: 

"Computational  Fluid  Dynamics  Is  making  an  Increasingly  major  Impact  In  aeronautical 
applications  and  on  the  aerodynamic  design  process.  The  rapid  progress  In  computer  capabil¬ 
ity,  the  general  availability  of  large  scale  computers,  and  the  parallel  achievements  In 
numerical  analysis,  algorithm  development  and  user  experience  assure  that  the  role  of  CFD  In 
aeronautics  will  continue  to  expand. 

“The  goal  of  the  Symposium  Is  to  provide  a  balanced.  If  not  exhaustive,  assessment  of  the 
status  of  CFD  In  aerodynamic  design  and  analysts.  The  emphasis  In  the  symposium  is  on  emerg¬ 
ing  applications  of  advanced  computational  techniques  to  complex  and  realistic 
configurations." 

The  Symposium  spanned  3-1/2  days  and  consisted  of  five  sessions  concentrating  on  four  major  topics: 
Session  I.  Grid  Generation 
Session  II.  Invlscld  Flow  I 
Session  III.  Invlscld  Flow  II 
Session  IV.  Vlscous-lnvlscld  Interactions 
Session  V.  Navler-Stokes 

As  suggested  by  the  titles  of  the  Individual  sessions,  the  speakers  addressed  numerous  detailed  aspects  of 
grid  generation  and  a  wide  range  of  numerical  methods  and  solutions.  The  titles  of  the  31  contributions 
are  listed  In  Appendix  A.  Most  of  the  papers  In  Session  II  were  concerned  with  potential  flows,  while 
Session  III  was  devoted  to  methods  for  solving  the  Euler  equations.  Both  Integral  ai  I  finite-difference 
boundary-layer  formulations  were  presented  In  Session  IV.  Most  of  Session  V  was  devoted  to  three- 
dimensional  solutions  of  the  Reynolds-averaged  Navler-Stokes  equations  with  turbulence  modeling  (abbre¬ 
viated  throughout  the  Symposium  as  simply  Navler-Stokes).  In  addition,  a  discussion  period  at  the  end  of 
the  Symposium  elicited  a  wide  range  of  Informal  comments  by  many  of  the  attendees.  The  transcript  of  this 
discussion  and  the  regular  papers  are  published  In  AGARO  Conference  Proceedings  CP-412.  The  members  of 
the  Program  Committee  for  the  Symposium  are  given  In  Appendix  B  of  this  report. 

In  keeping  with  the  stated  theme  of  the  Symposium  and  with  the  general  trends  In  the  aeronautical 
community  today,  most  of  the  papers  dealt  with  complex  configurations,  such  as  wing-body-store  combina¬ 
tions,  control  surfaces  within  jet  nozzles,  complete  fixed-wing  aircraft,  missiles,  and  helicopter  rotors 
and  fuselages.  In  general,  basic  fluid  dynamic  phenomena,  such  as  turbulence  modeling  or  direct  simula¬ 
tion  of  turbulence,  were  not  considered;  however,  fundamental  aspects  of  vortical  flows  and  tral 1 Ing-edge 
separation  were  treated  by  several  authors.  Also,  research  on  basic  numerical  analysis  and  new  algorithm 
development  was  less  emphasized  than  the  applications  of  existing/methods  and/or  their  extensions  to  three 
dimensions. 

It  may  be  mentioned  that  the  large  majority  of  the  papers  were  concerned  with  finite-difference  and 
finite-volume  methods,  apparently  reflecting  the  contemporary  tendencies  In  most  aeronautical  companies 
and  research  laboratories.  Also,  It  Is  noteworthy  that,  while  most  of  the  papers  emphasized  three- 
dimensional  (3-D)  flows,  only  three  papers  considered  unsteady  effects,  and  In  one  of  these  the  unsteadi¬ 
ness  In  the  solution  arose  as  an  unexpected  complication  In  a  problem  with  steady  boundary  conditions. 
However,  the  program  committee  had  not  solicited  papers  on  unsteady  flows,  as  these  had  been  addressed  In 
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two  previous  AGARD  meetings  [Refs.  1,2)  and  in  the  associated  reviews  [Refs.  3,4).  Finally,  hypersonic 
aerodynamics  was  discussed  by  only  two  speakers;  this  subject  is  the  basis  for  an  AGARD  Symposium  in 
April,  1987. 

2.  SYNOPSIS  OF  THE  PAPERS 

2.1.  Grid  Generation 

The  papers  In  Session  I  emphasized  the  great  challenges  that  grid  generation  poses  in  successfully 
applying  finite-difference  and  finite-volume  methods  to  meaningful,  practical  aeronautical  problems.  Even 
defining  the  surface  geometry  accurately  for  complex  3-D  bodies  becomes  a  time-  and  CPU-consuming  chore. 
Beyond  that,  creating  meshes  which  conform  to  3-D  body  surfaces  and  which  distribute  node  points  smoothly 
and  efficiently  throughout  flow  fields  Is  also  extremely  challenging  and  laborious,  especially  if  shock 
waves,  contact  discontinuities,  and/or  concentrated  vortices  are  present.  The  use  of  Interactive  computer 
graphics  is  indlspenslble,  and  automation  techniques  are  helpful,  but  considerable  human  Intervention  and 
ingenuity  Is  required  In  all  but  the  simplest  cases.  To  paraphrase  the  authors  in  this  session,  grid  gen¬ 
eration  for  complex  configurations  can  become  more  difficult  and  more  expensive  than  obtaining  the  numeri¬ 
cal  solution  itself.  And,  while  a  good  grid  will  not  guarantee  a  good  solution,  a  bad  grid  will  almost 
certainly  produce  a  bad  solution. 

PAPER  1.  OSKAM  and  HUIZING  described  a  2-D  zonal  grid  generation  method  that  Is  particularly  suitable 
for  multielement  airfoils.  Each  element  of  the  airfoil  is  embedded  in  one  or  more  separate  zones.  Within 
each  zone,  the  conformal  transformation  between  physical  and  computational  space  is  computed  by  a  varia¬ 
tional  formulation  that  minimizes  the  deviation  from  the  desired  grid  properties,  e.g.  cell  areas  and 
orthogonality,  while  allowing  points  to  move  along  some  of  the  zonal  boundaries.  This  freedom  of  movement 
facilitates  matching  the  grid  interfaces  between  zones,  and  local  surgery  can  be  done  to  relocate  topolog¬ 
ical  singularity  points  that  may  occur  at  zonal  corners.  An  example  was  given  of  a  grid  that  would  be 
suitable  for  Navler-Stokes  computations  for  an  alrfoll-f lap-slat  combination. 

PAPER  2.  WEATHERILL,  SHAW,  FORSEY,  and  ROSE  presented  a  zonal  technique  for  3-D  bodies  of  complicated 
geometry  such  as  complete  aircraft,  in  which  the  flow  domain  In  physical  space  is  subdivided  into  many 
r^nover lapping  "blocks."  A  feature  of  this  method  is  that,  in  principle,  the  grid  structure  and/or  topol¬ 
ogy  in  each  block  may  be  different  and  relatively  arbitrary,  according  to  what  is  most  appropriate  for 
that  zone.  However,  grid  lines  are  forced  to  match  and  pass  smoothly  across  the  interfaces.  The  grids 
are  generated  simultaneously  in  each  block  by  solving  sets  of  inhomogeneous  elliptic  partial  differential 
equations.  The  authors  candidly  discussed  some  of  the  problems  that  they  encountered,  such  as  grid  lines 
crossing  over  neighboring  lines,  the  difficulties  In  clustering  grid  points  arbitrarily,  and  slow 
convergence,  and  they  described  how  these  problems  were  ameliorated.  Results  were  shown  for  a  variety  of 
aircraft  and  missile  configurations  which  used  from  9  to  240  blocks.  The  computer  memory  required  to 
generate  a  complete  grid  system  is  about  half  that  needed  by  the  authors'  flow  solver  for  the  Euler  equa¬ 
tions,  whereas  the  CPU  time  is  generally  an  order  of  magnitude  less  to  generate  the  grid  than  it  is  to 
compute  the  solution. 

PAPER  3.  FRITZ  also  used  the  above  multiblock  concept  and  elliptic  partial  differential  equation 
(PDE)  solvers  in  his  3-D  grid  generation  method  for  complete  aircraft,  but  there  are  several  differences. 
His  grid  lines  are  continuous  across  zonal  boundaries,  but  the  slopes  are  not  matched.  Grid  points  are 
generated  in  each  block  In  steps;  first  along  the  perimeters,  then  on  each  surface,  and  finally  within 
each  block  separately.  However,  the  surface  grids  of  all  of  the  blocks  are  patched  together  before  the 
complete  grids  are  generated  inside  the  individual  blocks.  The  data  management  strategy  allows  large 
grids  of  several  million  grid  points  to  be  generated  on  computers  with  modest  main  memories.  However,  no 
figures  for  CPU  times  or  input/output  times  were  given. 

PAPER  4.  EDWARDS  described  both  a  procedure  for  defining  the  geometry  of  a  fighter  aircraft  and  a 
new,  efficient  grid-generation  method.  The  former  uses  an  existing  computer-aided  desiqn/computer-alded 
manufacturing  system  (CAD/CAM),  and  an  interactive  graphics  workstation  to  manipulate  the  input  data  base 
that  defines  seme  50  components  of  an  F-16  aircraft.  It  also  constructs  the  appropriate  surfaces  for  CFD 
modeling  of  the  wing,  fuselage,  canopy,  and  faired-over  inlet.  The  grid  systems  in  the  flow  field  are 
generated  by  applying  "parabolic"  difference  operators  in  two  of  the  three  directions,  to  the  elliptic 
PDEs  (without  the  source  terms)  used  by  the  previous  authors.  This  approximation  allows  a  noniterative, 
marching  solution  procedure  to  be  used  and  thereby  saves  considerable  CPU  time  and  in-core  memory,  while 
retaining  good  control  over  the  cell  areas,  clustering,  and  skewness.  In  the  example  given,  the  parabolic 
scheme  was  used  to  construct  a  relatively  coarse  grid  over  the  whole  flow  field,  and  finer  inner  grids 
were  interpolated  to  the  desired  distribution  on  the  body. 

2.2.  Invlscld  Flow  1 


This  session  was  mostly  devoted  to  linear  and  nonlinear  methods  for  potential  flows.  These  methods 
are  particularly  useful  In  the  aeronautical  design  process,  either  because  of  their  computational  effi¬ 
ciency,  or  for  their  ability  to  treat  complex  configurations  easily,  or  both. 

PAPER  5.  SMITH  and  WOODWARD  assessed  the  capabilities  and  deficiencies  of  three  established  panel 
methods  that  are  often  used  for  linear  calculations  of  complex  geometries  in  supersonic  flow.  They  com¬ 
pared  panel  solutions  for  relatively  simple  v  gs,  bodies,  and  wing-body  combinations  with  the  results  of 
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experiments,  higher-order  panel  methods,  and  Euler  calculations.  From  this  comparison,  they  Identified 
the  following  areas  where  problems  are  most  likely  to  arise: 

1.  Spurious  reflections  and  other  discrepancies  at  wing  tips, 

2.  Rounded  leading  edges  and  other  regions  of  high  surface  slopes, 

3.  The  neighborhood  of  wing-body  junctions. 

Rounded  leading  edges  that  produce  detached  bow  shock  waves  and  source  discontinuities  at  fuselage  panel 
edges  were  Indicated  to  be  the  most  troublesome. 

PAPER  6.  LE,  MORCHOISNE,  and  RYAN  described  a  panel  code  under  development  that  Incorporates  some  new 
numerical  techniques  to  speed  up  the  calculations.  They  have  developed  a  fast  Iterative  technique  using 
the  method  of  steepest  gradients,  and  a  special  way  of  introducing  extra  collocation  points  in  regions  of 
high  curvature.  This  method  was  shown  to  be  significantly  faster  than  a  direct  integration  technique, 
with  the  gain  increasing  as  the  number  of  panels  Increased.  Some  preliminary  results  were  shown  for  a 
sphere,  wing,  and  helicopter  and  transport  fuselages. 

PAPER  7.  STRAWN  and  TUNG  discussed  the  adaptation  of  an  unsteady,  3-D  full-potential  code  to  compute 
the  transonic  flow  near  the  tips  of  helicopter  rotor  blades.  For  this  class  of  problems,  the  vortical 
wake  Is  very  important,  but  calculating  its  structure  and  position  explicitly  or  directly  Is  beyond  the 
current  state  of  the  art  In  CFD.  Therefore,  two  different  methods  were  developed  to  couple  approximate 
rotor  wake  models  with  the  finite-difference  code.  The  first  formulation  Iterates  between  (1)  the  local 
blade  lift  which  Is  calculated  by  the  full-potential  code  for  a  given  angle  of  attack  distribution,  and 
(2)  the  blade- surface  Inflow,  or  effective  angle  of  attack,  which  is  obtained  from  an  approximate.  Inte¬ 
gral  wake  method.  The  Integral  code  also  predicts  the  instantaneous  blade  motion.  The  second  wake¬ 
modeling  scheme  allows  concentrated  vortices  of  predetermined  strength  to  be  Introduced  into  the  interior 
of  the  computational  domain  by  means  of  a  "split-potential"  formulation  of  the  finite-difference  equa¬ 
tions.  Although  this  work  Is  not  complete,  encouraging  agreement  with  experiments  was  obtained  for  both 
methods. 

PAPER  8.  KAFYEKE  presented  a  3-D  transonic  small-disturbance  code  that  uses  grid  embedding  to  predict 
the  aerodynamic  interference  of  wing-body-pylon-store  configurations.  The  mean-surface  approximation 
inherent  in  the  small -disturbance  formulation  greatly  simplifies  the  grid-generation  problem  by  allowing 
cartesian  meshes  to  be  used  everywhere,  without  conforming  to  the  actual  surfaces.  Furthermore,  fine 
grids  around  each  component  are  embedded  within  an  overall  coarse  grid,  and  the  solution  procedure  cycles 
between  the  coarse  and  fine  grid  regions  to  accelerate  the  convergence.  Citing  the  good  agreement 
obtained  with  experiments  and  the  modest  computational  requirements  of  the  code  when  using  approximately 
200,000  grid  points,  the  author  downplayed  the  need  for  supercomputers  In  the  aeronautical  design  process. 

PAPER  9.  VAN  DER  VOOREN,  VAN  DER  WEES,  and  MEELKER  described  the  computational  aerodynamics  Inte¬ 
grated  system  called  MATRICS  (Multi-component  Aircraft  Transonic  Inviscld  Computational  System)  that  Is 
under  development  at  the  National  Aerospace  Laboratory.  This  system  currently  solves  the  full -potential 
equation  using  a  multigrid  method  enhanced  at  National  Aerospace  Laboratory  NLR,  with  provisions  for 
approximately  modeling  vortex  sheets,  propeller  slipstreams,  and  jet  exhaust  plumes.  Extensions  are 
underway  to  the  Euler  equations  in  subdomains,  and  eventually  to  the  Navier-Stokes  equation.  MATRICS  is 
designed  with  both  scalar  and  vector  computers  in  mind,  but  all  the  data  are  stored  in  the  main  memory 
without  I/O  (Input/output)  transfer.  In  anticipation  of  the  large  memories  that  future  supercomputers  will 
feature.  In  addition  to  the  systems  aspects,  the  authors  reported  on  effects  of  boundary  conditions, 
free-stream  "consistency,"  artificial  viscosity,  and  the  size  of  the  computational  domain  on  the  computed 
results  for  a  representatl ve  transonic  wing  with  and  without  a  cylindrical  body  and/or  propeller.  They 
noted  that  the  calculation  of  drag  entails  more  stringent  requirements  for  both  grids  and  convergence 
tolerences  than  does  lift. 

PAPER  10.  PETRIE  and  SINCLAIR  gave  a  progress  report  on  their  development  of  a  nonlinear  field  panel, 
or  integral,  method  for  solving  the  ful 1 -potential  and  Euler  equations  for  complex  configurations.  In 
this  approach,  field  computational  grids  are  only  needed  In  relatively  small  regions,  and  they  need  not  be 
body-conforming.  For  compressible  flows,  however,  nonlinear  volume  Integrals  arise  that  must  be  evaluated 
Iteratively.  No  details  of  the  solution  technique  nor  of  the  computational  requirements  and  efficiency 
were  given. 

PAPER  11.  MARCHBANK  presented  an  overview  of  the  use  of  CFD  in  military  aircraft  design,  Including 
several  examples  of  fighter  aircraft  In  supersonic  flow.  He  stressed  three  Important  requirements  for 
effective  use  of  CFD: 

1.  Speed  -  fast  interaction  with  the  aircraft  geometry  definition, 

2.  Utility  -  easy-to-use  CFD  procedures,  oriented  to  the  engineer, 

3.  Credibility  -  adequate  accuracy  of  the  aerodynamic  results. 

Panel  methods  more  nearly  satisfy  the  first  requirement,  both  with  regard  to  solution  time  and  to  geometry 
Interface,  whereas  Euler  calculations  were  found  to  be  significantly  more  accurate  for  the  supersonic 
flows  of  Interest.  Examples  were  given  of  drag  optimization,  wing  deslgi.  with  drag  constraints, 
pi tch log -moment  optimization,  and  forebody  and  canopy  design.  These  results  were  shown  to  lead  to  signif¬ 
icant  improvements  In  supersonic  aircraft  performance  without  compromising  transonic  characteristics, 
within  short  project  timescales. 

PAPER  12.  W<  \DL0W  and  DAVIS  began  the  sequence  of  papers  on  Euler  methods  that  continued  through 
Session  III.  Their  method  is  a  space-marching,  finite-volume  Implementation  of  a  Godunov-type  scheme  In 
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supersonic  flow.  It  should  be  explained  that  more  conventional  schemes  assume  a  smooth  (or  mostly 
smooth),  flow  field,  whereas  the  Godunov  method  considers  the  flow  to  consist  of  a  series  of  piece-wise 
constant  states  with  discontinuities  occurring  midway  between  mesh  points  (Ref.  5j.  Analytic  properties 
of  2-0  shock  waves  and  expansion  fans  are  used  to  advance  the  solution  across  cell  faces.  In  this  paper, 
the  addition  of  a  central-difference  predictor  step  and  slope  limiters  gives  second-order  accuracy.  This 
approach  Is  more  robust  without  explicitly  adding  artificial  viscosity  than  are  many  conventional  finite- 
difference  schemes,  although  more  computational  work  is  required  per  grid  point.  Reasonable  results  were 
reported  for  a  variety  of  finned  tactical  missile  configurations. 

2-3.  Invlscld  Flow  II 

Euler  codes  are  emerging  rapidly  from  the  research  stage,  and  as  demonstrated  in  this  session,  they 
are  finding  increasing  applications  to  3-D  configurations  with  strong  nonlinear  features.  A  variety  of 
methods  with  a  wide  range  of  computer  requirements  were  presented  for  solving  the  Euler  equations. 

PAPER  13.  BRE01F,  CHATTOT,  KOECK,  and  WERLE  used  an  explicit,  multigrid  finite-volume  method  to  cal¬ 
culate  the  transonic  and  supersonic  flow  in  an  axlsymnetric  nozzle  with  control  surfaces  protruding  into 
the  stream  just  upstream  of  the  nozzle  exit.  A  relatively  coarse  grid  system  of  approximately  60,000 
nodes  In  two  blocks  was  used,  requiring  on  the  order  of  10  min.  CPU  time  on  a  Cray  IS  computer.  The 
results  were  compared  with  linear  theory  and  experiments;  It  was  found  that  the  thrust  deflection  caused 
by  the  control  surface  was  accurately  predicted,  but  the  loss  in  thrust  which  was  due  to  the  interaction 
was  underestimated. 

PAPER  14.  LEICHER  employed  a  different  explicit,  multiblock,  multigrid  finite-volume  method  to  calcu¬ 
late  transonic  flow  over  a  wing  and  subsonic  flow  in  the  diffusers  of  a  turbine  and  a  wind  tunnel.  A 
wing -propeller  combination  was  also  simulated  using  an  actuator-disk  model  at  the  plane  of  the  propel¬ 
ler.  Extensive  grid-refinement  and  mesh-spacing  studies  were  carried  out  for  the  wing,  with  the  number  of 
grid  points  varying  from  approximately  20,000  to  1,200,000.  From  the  Information  given  in  the  paper,  the 
latter  appeared  to  require  over  100  hours  of  CPU  time  on  an  IBM  3083  computer,  but  an  intermediate  grid  of 
about  150,000  points  seemed  to  give  satisfactory  results  In  much  less  time.  The  results  for  the  Kaplan 
turbine  diffuser  showed  evidence  of  bimodal  steady-state  solutions. 

PAPER  15.  MURMAN  and  RI2ZI  provided  one  of  the  more  fundamental  and  provocative  presentations  of  the 
symposium.  Their  paper  concerned  the  capabilities  of  Euler  codes  to  generate  and  simulate  leading-edge 
vortices  on  delta-wing  configurations  at  high  Reynolds  numbers.  Several  examples  were  given  of  delta 
wings  with  sharp  leading  edges,  in  both  subsonic  and  supersonic  flow.  Concentrated  vortical  structures 
emerged  in  the  solutions  with  approximately  the  correct  total  pressure  loss  In  the  centers  of  the 
structures.  From  their  own  and  other  studies,  they  found  this  total  pressure  loss  to  be  approximately 
Independent  of  numerical  parameters,  such  as  grid  spacing  and  the  level  of  the  artificial  dissipation. 
Also,  the  calculated  total  pressure  loss  was  relatively  Independent  of  whether  or  not  physical 
viscous-stress  terms  were  Included  in  the  equations.  On  the  other  hand,  this  loss  was  found  to  be  sensi¬ 
tive  to  flow  and  geometry  parameters,  as  It  should  be.  Professor  Murman  explained  his  argument  that 
“...any  mechanism,  whether  real  or  artificial  (i.e.,  numerical  -  author),  which  gives  the  vortex  sheet 
some  thickness  will  lead  to  a  total  pressure  loss."  He  went  on  to  discuss  the  structure  of  vortex  cores, 
some  Implications  of  total  pressure  losses  with  respect  to  vortex  breakdown,  and  solutions  that  were  glob¬ 
ally  stationary  but  locally  unsteady.  Overall,  this  paper  Identified  phenomena  which  will  manifest  them¬ 
selves  In  future  calculations  and  which  will  surely  be  the  subject  of  further  controversy  and 
investigation. 

PAPER  16.  PERRIER,  whose  group  at  Avions  Marcel  Dassaul t-Breguet  has  pioneered  the  use  of  finite- 
element  methods  for  aeronautical  applications,  gave  a  brief  review  and  comparison  of  f inite-volume  and 
finite-element  methods  for  solving  the  Euler  equations.  He  expressed  the  view  that  Navler-Stokes  solvers 
are  too  expensive  for  routine  industrial  use,  but  he  stressed  that  potential  methods  fall  to  predict  many 
nonlinear  aerodynamic  phenomena  accurately  enough.  He  showed  numerous  examples  of  solutions  for  complex 
flows  and/or  configurations,  and  he  discussed  the  role  of  Euler  codes  as  approximations  to  Navier-Stokes 
codes  and  as  possible  aids  in  understanding  turbulence. 

PAPER  17.  EBERLE  and  MISEGADES  described  a  new  finite-volume  Euler  code  that  is  claimed  to  be  third- 
order  accurate  and  which  can  be  run  In  either  an  explicit  or  Implicit  mode.  Their  method  also  uses 
Godunov  concepts  (see  Paper  12,  above).  A  relatively  simple  grid-generation  scheme  with  grid  embedding 
allowed  them  to  generate  several  Impressive  solutions,  using  up  to  520,000  cells,  for  a  complete  fighter 
aircraft  In  transonic  and  supersonic  flow  (see  Fig.  1.  This  figure  was  In  color  In  the  origlna-'  paper). 

A  limited  grid-ref Inement  study  Indicated  that  grids  of  this  size  were  required  for  reliable  results.  The 
CPU  times  for  the  computations,  which  were  done  on  Cray  X-MP  machines,  were  not  quoted,  but  based  on  the 
number  of  grid  points  and  the  properties  of  the  code,  one  might  Infer  that  they  were  several  hours. 

PAPER  18.  KARMAN,  STE IN8RENNER,  and  KISIELEWSKI  also  presented  a  demonstration  Euler  calculation  for 
a  complete  fighter  aircraft,  although  their  solution  was  not  not  fully  converged  at  the  time  of  the  Sympo¬ 
sium.  They  started  with  a  multiblock,  explicit  finite-volume  predictor-corrector  formulation,  and 
employed  20  blocks  with  a  total  of  about  530,000  grid  points.  After  1500  Iterations,  they  switched  to  an 
implicit  flow-solver  with  a  faster  rate  of  convergence,  but  this  cheme  was  not  described.  The  results 
presented  after  3100  Iterations  appeared  to  be  evolving  toward  agreement  with  experimental  data.  In  a 
private  cownunlcatlon,  the  authors  have  since  confirmed  this  trend  after  another  1000  Iterations,  although 
the  experimental  multiple-shock  structure  on  the  wing  never  materialized  In  the  solution;  this  behavior  s 
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Fig.  1.  Computed  surface  pressure  distributions  on  a  fighter  aircraft;  Paper  17. 

attributed  to  inadequate  grid  resolution.  The  code  was  run  on  Cray  X-MP  and  Cray  2  machines  for  a  total 
of  some  30-35  hours,  of  which  approximately  30*  was  I/O  time  to  the  SDD  (solid-state  device)  external 
memory  device. 

2.4  Viscous-Inviscid  Interactions 

As  suggested  by  the  title,  the  papers  in  Session  IV  addressed  the  relatively  inexpensive  introduction 
of  viscous  effects  by  means  of  boundary- layer  concepts  coupled  with  invfscld  flow  solvers.  Concerning 
this  approach,  a  brief  explanation  of  the  nomenclature  used  in  this  report  may  be  In  order.  The  coupling 
is  called  ‘'weak''  If  the  inviscld  pressure  distribution  Is  Impressed  directly  on  the  boundary  layer  and  the 
resultant  boundary  layer  weakly  perturbs  the  Inviscld  solution.  At  least  in  2-0  steady  flows,  this  proce¬ 
dure  becomes  singular  in  the  neighborhood  of  flow  reversal.  "Strong"  coupling  means  that  the  boundary- 
layer  calculation  Is  performed  in  an  Inverse  mode  and  the  computed  pressure  distribution  Is  intimately  and 
strongly  coupled  with  the  inviscld  flow-solver.  Many  of  the  special  techniques  that  have  been  introduced 
to  accelerate  convergence  or  to  Improve  stability  tend  to  blur  this  formal  distinction  between  weak  and 
strong  coupling,  but  the  “strong"  terminology  Is  retained  herein  when  their  Intent  is  to  achieve  strong, 
coupling  results. 

PAPER  19.  REIS  and  THOMPSON  compared  weakly-coupled  calculations  by  two  different  numerical  methods 
with  measurements  of  2-0  tralllng-edge  separation.  Their  study  Included  a  laudable,  careful  assessment  of 
numerical  errors  and  uncertainties.  However,  they  found  that  while  the  numerical  errors  were  different, 
the  discrepancies  with  the  experiment  were  about  the  same;  therefore,  the  turbulence  model  was  blamed. 

Not  surprisingly,  computed  surface  pressures  (lift)  were  more  satisfactory  than  skin  friction  and  wake 
profiles  (drao). 

PAPER  20.  SCHMAT2  and  H1R5CHEL  presented  two  zonal  schemes  for  coupling  the  Euler  equations  with  both 
the  boundary- layer  and  Navler-Stokes  equations,  depending  on  the  region  of  the  flow  field  and  the  strength 
of  the  Interaction.  In  the  "alternating  (weak)  coupling"  procedure  of  cyclic  Iteration  through  three 
distinct  zones,  different  numerical  schemes  were  used  for  each  equation  set,  with  the  grids  overlapping 
between  the  Euler  and  Navler-Stokes  zones.  Good  results  were  shown  for  airfoils  without  separation.  In 
the  “close  coupling"  method,  weak  Interaction  of  the  boundary  layer  was  retained  on  the  front  of  the  air¬ 
foil,  but  the  Navler-Stokes  grid  near  and  behind  the  body  was  embedded  In  an  outer  Euler  grid.  The  Euler 
and  Navler-Stokes  zones  were  solved  simultaneously  with  a  common  implicit  relaxation  algorithm.  In  2-0,  a 
speedup  of  about  a  factor  two  was  realized  compared  to  a  full  Navler-Stokes  solution;  the  gain  in  3-0 
problems  remains  to  be  determined. 


PAPER  21  was  withdrawn,  regrettably. 

PAPER  22.  GULCAT  studied  a  3-0  unsteady  boundary  layer  on  a  body  placed  Impulsively  into  motion.  He 
solved  the  Incompressible  laminar  boundary- layer  equations  in  physical  space,  eschewing  similarity  trans¬ 
formations.  He  employed  a  streamwise-  and  time-marching  scheme  that  Is  Implicit  In  the  normal  direction; 
the  method  seems  to  follow  the  spirit  of  Blottner's  method  (Ref.  6|  for  2-0  steady  boundary  layers.  The 
resultant  code  Is  relatively  simple,  computationally  efficient,  and  stable  for  nonsingular  boundary- layer 
flow.  He  obtained  good  results  for  the  flow-reversal  time,  pressure  distribution,  and  pressure  distribu¬ 
tion  for  the  Impulsive  flow  past  a  circular  cylinder.  Solutions  on  the  plane  of  symmetry  of  an  oblate 
spheroid  were  also  given  up  to  the  time  of  the  development  of  flow  reversal  at  the  rear  stagnation  point. 

PAPER  23.  CLER  presented  a  number  of  3-D  solutions  obtained  on  helicopter  fuselages,  air  intakes,  and 
hub  fairings.  This  effort  was  largely  motivated  by  the  need  to  estimate  the  drag  of  bluff  bodies  with 
large  separated  base  flows.  Most  of  the  results  were  obtained  with  a  panel  code  combined,  but  not 
directly  coupled  with,  a  3-0  Integral  boundary- layer  code  that  located  the  lines  of  separation.  While 
obviously  not  reproducing  all  the  flow  features,  these  solutions  often  help  Identify  troublesome  regions 
and  suggest  modifications  to  the  body  geometry.  The  crucial  need  for  an  adequate  vortex-wake  model  was 
noted,  and  preliminary  results  were  shown  that  Had  been  obtained  with  the  0NERA  discrete-vortex  method  for 
a  rotating  blade  In  forward  flight.  However,  this  technique  was  found  to  be  difficult  to  combine  with  the 
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fuselage  panel  and  boundary- layer  programs,  and  much  more  development  will  be  required  to  meet  this 
objective. 

PAPER  24.  VAN  DAlSEM  and  STE6ER  described  a  simple  and  efficient  algorithm  for  solving  the  unsteady, 
3-0  boundary- layer  equations  In  either  a  time-accurate  or  relaxation  mode.  This  approach  is  claimed  to  be 
more  flexible  and  easy  to  apply  than  space-marching  procedures.  Their  code  switches  from  a  direct  (weak) 
to  an  Inverse  (strong)  mode  near  and  within  reverse-flow  regions,  which  permits  separation  to  be  computed 
readily.  Two  examples  of  good  results  were  shown  for  prescribed  pressure  distributions.  They  also  exam¬ 
ined  the  use  of  this  boundary- layer  algorithm  to  speed  up  the  convergence  of  an  existing  Navler-Stokes 
flow  solver.  The  combination  results  In  a  so-called  "Fortified  Navler-Stokes"  scheme  that  looks  very 
promising.  Fine-grid  boundary- layer  results  near  the  wall  are  used  as  forcing  functions  In  the  thin-layer 
Navler-Stokes  equations,  which  are  solved  only  on  a  coarse  grid.  Two  examples  were  given  in  which  the 
accuracy  of  fine-grid  Navler-Stokes  solutions  was  obtained  20  times  faster  with  the  Fortified  Navler- 
Stokes  scheme. 

PAPER  25.  LAZAREFF  and  LE  BALLEUR  described  their  Multi-Zonal -Marching  (MZM)  method  for  turbulent, 

3-0  boundary  layers  with  vlscous-lnviscld  Interaction.  The  MZM  method  Is  a  procedure  for  extending 
Le  Balleur's  well-established  and  successful  2-0  Integral  "semi- Inverse"  (essentially  strong  coupling) 
techniques  to  three  dimensions.  Careful  attention  Is  given  to  the  characteristic  cones  of  Influence  in 
each  zone,  and  multiple  sweeps  can  correctly  couple  different  flow  domains  with  a  wide  range  of  flow 
directions.  The  computations  of  swept  wings,  for  example,  start  at  a  leading-edge  stagnation  point,  from 
which  multiple  sweeps  are  made  along  the  leading  edge,  but  within  a  narrow  streamwlse  zone,  before  the 
final  streamwlse  and  spanwise  sweeps  progress  over  the  upper  and  lower  surfaces.  Limited  results  were 
presented  of  a  wing  calculation  with  strong  vlscous-lnviscld  coupling,  but  with  no  flow  separation.  How¬ 
ever,  most  of  the  results  were  obtained  in  the  direct  boundary- layer  mode  for  ellipsoids  and  a  reentry- 
type  lifting  body  at  angle  of  attack.  The  MZM  strategy  allowed  the  calculations  to  be  continued,  with  no 
vlscous-lnviscld  coupling,  over  the  whole  body,  even  though  crossflow-separation  behavior  developed  on 
these  bodies. 

PAPER  26.  FERMIN  concluded  this  session  with  an  outline  of  two  RAE  viscous-interaction  wing  codes  and 
a  description  of  several  swept -wing  applications.  Transonic  small-disturbance  and  full  potential  codes 
have  been  coupled  with  a  3-0  integral  boundary- layer  method  in  the  direct  mode.  The  former  code  was  used 
to  design  the  new  RAE  M2155  low-aspect-ratio  research  wing  that  comprises  a  major  code-validation  exercise 
In  progress,  and  wind  tunnel  results  from  the  Initial  tests  of  this  wing  were  presented.  Fairly  detailed 
comparisons  were  made  with  this  experiment  and  with  one  of  a  transport-wing  and  body  combination.  Compu¬ 
tations  done  with  the  Viscous  Full  Potential  code,  with  direct  (weak)  coupling,  were  In  excellent  agree¬ 
ment  with  the  measurements  In  the  subsonic  case.  However,  transonic  conditions  with  strong 
vlscous-lnviscld  interaction  were  a  major  challenge.  It  was  difficult  to  obtain  converged  solutions,  and 
the  computed  pressure  distributions  were  deficient  In  local  regions,  particularly  when  the  boundary  layer 
was  close  to  separating.  Suggestions  for  Improvements  Included  a  different  coupling  scheme,  the  use  cf  an 
Inverse  (strong)  method,  and  using  Euler  codes  for  the  Invlscld  flow. 

2.5  Navier  Stokes 

The  increasing  availability  of  supercomputers  and  the  progress  In  algorithm  and  software  development 
have  brought  3-0  Navler-Stokes  calculations  of  aeronautical  configurations  close  to  the  realm  of  near-term 
reality.  The  papers  In  this  session  gave  valuable  insights  Into  the  general  state  of  the  art  today,  and 
Into  the  developing  trends  In  applied  computational  aerodynamics  for  viscous  problems. 

PAPER  27.  SHANG  outlined  the  most  common  approaches  used  In  the  aerospace  community  and  reviewed  a 
number  of  recent  accomplishments.  Including  simulations  of  flows  3»*ound  aircraft  wings,  fuselages,  after¬ 
bodies,  and  Inlets.  He  briefly  discussed  numerical  efficiency  ar j  accuracy,  boundary  conditions,  turbu¬ 
lence  models,  grid  generation,  data  structure  and  management,  and  post-processing  and  display  of 
results.  He  noted  that  the  grid-generation  phase  occupies  an  ever- Increasing  fraction  of  the  total 
elapsed  time  of  numerical  simulations,  and  he  lamented  the  vast  amount  of  information  that  is  discarded 
during  or  soon  after  a  typical  investigation.  The  growing  ana  lasting  importance  of  "Interdisciplinary 
computational  fluid  dynamics,"  In  which  Navler-Stokes  equations  are  coupled  with  the  governing  equations 
of  solid  mechanics,  chemistry,  combustion,  electromagnetics,  optics,  etc.,  was  also  stressed. 

PAPER  28.  WA1,  BLOW,  and  Y0SH1HARA  employed  the  Parabolized  Navler-Stokes  (PN$)  method.  In  which  the 
streamwlse-d If fusion  terms  In  the  equations  are  neglected,  to  calculate  the  supersonic  flow  over  a  generic 
fighter  aircraft  configuration  at  high-lift  conditions.  They  used  an  existing  PNS  code,  described  In  the 
next  paper,  with  local  regions  of  subsonic  flow  computed  by  an  unsteady  Navler-Stokes  code.  Oata  for 
comparison  were  very  sparse,  but  reasonable  agreement  of  the  numerical  and  experimental  results  was 
obtained.  The  calculation  time  for  approximately  4100  grid  points  In  the  crossflow  planes  and  about  1000 
streamwlse  stations  was  5.4  hours  on  an  ETA  205  computer  for  a  wing,  fuselage,  canard,  and  nacelle  combi¬ 
nation.  A  major  fraction  of  this  time  was  consumed  In  regenerating  the  grid  at  each  streamwlse  station. 

PAPER  29.  CHAUSSEE  discussed  the  NASA-Ames  PNS  code  and  several  applications  to  supersonic  and  hyper¬ 
sonic  configurations  with  Important  viscous  effects.  The  parabolized  approximation  assumes  that  the 
Invlscld  flow  is  supersonic  In  the  streamwlse  direction,  and  that  the  subsonic  flow  In  the  viscous  sub¬ 
layer  Is  always  positive  In  the  streamwlse  direction.  This  allows  the  solution  to  march  streamwlse,  giv¬ 
ing  substantial  savings  In  CPU  time,  while  an  Implicit  algorithm  Is  used  to  solve  the  flow  In  each  cross- 
flow  plane.  Crossflow  separation  Is  permitted,  but  streamwlse  separation  Is  not.  An  algebraic 


eddy-viscosity  model  Is  used  In  turbulent  cases.  Examples  Included  ogive-cylinders,  sphere-cones  with 
flaps,  finned  projectiles,  the  Space  Shuttle  orblter,  and  the  supersonic  flgbcer  which  is  the  subject  of 
Paper  28. 
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PAPER  30,  Part  I.  FLORES,  HOLST,  KAYNAK,  GUNDY,  and  THOMAS  described  the  development  of  a  zonal 
Euler/Navler-Stokes  code  called  Transonic  Navier-Stokes  (TNS).  The  method  uses  a  diagonalized  implicit 
algorithm  to  solve  the  unsteady  equations,  but  with  a  variable  time  step  for  faster  convergence.  They 
divided  the  flow  Into  four  zones  for  wings:  two  Inner  blocks  near  the  surface  where  the  thin-layer 
Navier-Stokes  equations  are  solved,  and  two  outer  Euler  blocks.  As  In  the  preceding  two  papers,  turbu¬ 
lence  Is  simulated  by  means  of  an  algebraic  eddy-viscosity  model.  Results  for  an  airfoil  in  a  solid-wall 
wind  tunnel  and  for  two  transonic  wings  were  given;  the  agreement  with  experiments  was  generally  good, 
although  some  differences  in  regions  of  shock-induced  separation  were  noted.  CPU  times  on  a  Cray  X-MP 
computer  were  on  the  order  of  one  hour  for  150,000  grid  points. 

PAPER  30,  Part  II.  CHADERJIAN  applied  the  aforementioned  TNS  code  to  a  swept  NACA  0012  wing  at  high 
angle  of  attack,  with  a  view  toward  demonstrating  the  robustness,  efficiency,  and  accuracy  of  the  code, 
and  assessing  the  characteristics  of  the  grid.  Major  improvements  near  grid  singularities  were  accom¬ 
plished  by  computing  the  transformation  metrics  consistently  to  preserve  a  uniform  free  stream.  Results 
were  obtained  up  through  maximum  lift  at  M  =  0.5  and  0.8,  although  data  were  not  available  for  compari¬ 
son.  The  author  suggested  Improvements  In  terms  of  additional  grid  refinement,  an  improved  numerical- 
dissipation  model,  and  a  turbulence  model  more  suited  to  shock-induced  separation. 

PAPER  31.  KORDULLA,  VOLLMERS,  and  DALLMANN  reported  on  their  simulation  of  transonic  laminar  flow 
past  a  hemisphere-cylinder  at  angle  of  attack,  and  they  analyzed  the  topology  of  the  separated  flow  in 
great  detail.  Their  time- accurate,  finite-volume  method  uses  a  variation  of  the  MacCormack  explicit- 
implicit  predictor-corrector  scheme,  and  up  to  about  230,000  grid  points  were  used  in  the  present  study; 
CPU  times  on  the  Cray  IS  computer  used  were  not  quoted.  The  fine-mesh  solution  of  a  nominally  steady  case 
exhibited  unsteadiness  in  the  separated  region,  and  the  question  of  whether  this  was  due  to  fluid  physics 
or  to  numerical  errors  was  not  resolved.  The  main  thrust  of  the  paper,  however,  was  the  analysis  of  the 
topology  of  the  computed  flow  field.  Guided  by  earlier  work  of  the  third  author  and  others  on  topological 
kinematics  and  bifurcation  theory,  considerable  detail  was  extracted  from  the  solution.  Impressive 
sketches  of  complex  vortical  structures  and  of  the  major  singular  stream  surfaces  were  presented  (see 
Figs.  2  and  3). 


Fig.  2.  Sketch  of  the  major  singular  points  In  the  symmetry  plane  and  on  the  wall  of  a  hemispnere 
cylinder  at  angle  of  attack  in  transonic  flow;  Paper  31. 


Fig.  3.  Sketch  of  the  major  singular  stream  surface  In  the  flow  field  of  a  hemisphere  cylinder; 

Paper  31. 


PAPER  32.  Finally,  SHANG  discussed  his  pioneering  computations  of  the  flow  field  of  the  complete 
XC-240  experimental  aircraft  at  a  Mach  number  of  5.95.  The  equations  In  weak  conservation- law  form  were 
solved  In  a  single  computational  block  using  MacCormack' s  explicit  unsplit  algorithm,  with  the  aid  of  a 
sophisticated  data  management  scheme  and  the  Solid  State  Device  (SSD)  of  the  Cray  X-MP.  The  solution  time 
for  475,200  grid  points  was  about  20  hours,  of  which  about  25*  was  consumed  by  I/O  transfer  to  the  SSD. 
Good  agreement  was  obtained  between  computed  and  measured  lift  and  drag,  and  on  the  whole,  the  pressure 


distributions  and  heat  transfer  were  satisfactory.  A  major  effort  to  display  the  results  with  the^ald  of 
3-D  color  graphics  produced  some  spectacular  pictures  of  the  surface  streamlines,  pressures,  and 
temperatures. 

3.  DISCUSSION  AND  EVALUATION 

The  goal  of  the  Symposium,  as  noted  In  the  Introduction,  was  to  provide  an  ”...  assessment  of  the 
status  of  CFO  in  aerodynamic  design  and  analysis,"  with  an  emphasis  on  "...  emerging  applications  of 
advanced  computational  techniques  to  complex  and  realistic  configurations.”  Indeed,  the  material  pre¬ 
sented  at  the  Symposium  was  generally  indicative  of  the  state  of  the  art  of  applied  computational  aerody¬ 
namics  In  the  NATO  countries,  in  terms  of  current  activities,  methods,  numerical  and  physical  modeling 
employed,  and  computers  In  service. 

However,  the  Symposium  had  definite  limits  in  scope  and  content.  The  Program  Committee  was  obliged 
to  leave  out  noteworthy  efforts  in  almost  every  country,  and  there  was  little  representation  from  Sweden 
and  none  from  Japan.  Also,  research  on  new  algorithm  development  and  turbulence  modeling  was  generally 
excluded,  apparently  by  Intent.  Beyond  that,  It  was  somewhat  surprising  that  the  work  of  Anthony  Jameson 
and  his  colleagues  on  the  FLO-xx  series  of  codes  and  their  many  permutations  was  hardly  represented. 

As  a  general  consent  on  Session  IV,  Viscous-Inviscid  Interactions,  weakly-coupled  3-D  methods  have 
been  developed  In  several  countries  over  the  past  few  years,  but  these  efforts  were  not  well  repre¬ 
sented.  It  Is  not  clear  whether  this  was  accidental,  or  whether  it  reflects  a  new  trend  toward  developing 
strong -coupl ing  methods  that  simply  has  not  yet  materialized. 

The  Symposium  did  make  clear  the  widespread  use  of  CFD  throughout  the  aeronautical  Industry,  and  It 
revealed  that  the  spectrum  of  methods  in  use  is  growing  rapidly  to  include  Euler  and  Navier-Stokes  codes 
applied  to  complex  configurations.  The  Symposium  was  also  very  timely,  because  the  emphasis  has  clearly 
shifted  from  2-D  problems  to  the  more  realistic  world  of  3-0.  And  as  discussed  in  Section  3.1.2  below, 
the  advent  of  this  geometrical  complexity  is  causing  the  time,  manpower,  and  computer  costs  of  grid  gener¬ 
ation  and  data  analysis  to  become  much  more  significant  relative  to  the  effort  previously  involved  In 
simply  generating  solutions. 


3. 1  Costs  and  Capabilities 

Today  the  aerodynamic  1st  has  a  wide  range  of  computational  methods  available,  at  a  correspondingly 
wide  range  in  manpower  and  computer  costs.  As  mentioned  in  the  Introduction,  finite-difference  or 
finite-volume  methods  remain  much  more  popular  for  nonlinear  problems  than  for  field-panel  or 
finite-element  methods.  M.  Pierre  Perrier  argued  persuasively  for  the  latter,  but  no  direct  quantitative 
comparisons  were  made,  nor  could  be  made  from  the  material  presented  at  the  Symposium.  For  inviscid 
flows,  3-D  potential  codes  continue  to  be  developed  and  improved,  but  this  approach  is  now  relatively 
mature  and  the  development  of  Euler  methods  has  become  more  fashionable.  However,  the  treatment  of  vis¬ 
cous  effects  has  not  kept  pace,  at  least  in  3-0.  Vlscous/lnviscld  coupling  methods  are  potentially  very 
cost-effective  in  CPU  time,  but  much  remains  to  be  done  to  make  them  work  properly  for  flows  with  strong 
shock  waves.  Pioneering  showcase  Navier-Stokes  solutions  are  emerging  rapidly;  but  as  noted  below,  the 
CPU  times  remain  very  long,  the  convergence  is  slow,  and  questions  about  accuracy,  the  effects  of  grid 
refinement,  and  turbulence  modeling  have  barely  been  asked,  let  alone  answered. 


3.1.1  Computational  Requirements.  The  different  levels  of  approximations  to  the  governing  equations 
of  fluid  mechanics,  and  the  associated  difficulties  In  solving  them  numerically,  are  well  known  and  need 
not  be  reiterated  here.  However,  it  is  Instructive  to  examine  briefly  the  dominant  factors  that  determine 
the  computational  requirements  for  the  range  of  methods  that  were  presented  at  the  Symposium.  Of  course, 
many  factors  such  as  the  degree  of  maturity  of  different  algorithms  and  the  Impact  of  new  computer  archi¬ 
tectures  are  difficult  to  quantify;  but  the  approximate  CPU  time  required  for  most  of  the  methods  can  be 
crudely  estimated  with  the  aid  of  the  following  formula: 


CPU  =  A  .  Mg  .  Ng  .  Nj/FlOPS 


where: 


A 


ng  * 

m 

*1  = 

flops  = 


"numerical  inefficiency" 
number  of  operations  per 
number  of  grid  points  or 
1  for  finite-difference 
number  of  iterations,  or 
number  of  floating-point 


factor 

grid  point  per  time  step 
panels 

,  »  2  for  panel  methods 

number  of  time  steps  for  an  unsteady  calculation 
arithmetic  operations  per  second 


(1) 


The  inefficiency  factor.  A,  Is  Introduced  here  to  emphasize  that  few  codes  take  full  advantage  of  the 
computer  being  used;  in  practice  this  factor  Is  a  function  of  the  programing  efficiency,  the  degree  of 
vectorization,  the  coupling  between  the  grid  and  the  solution  algorithm,  the  data-management  strategy,  the 
user  experience,  etc.  Ideally,  Its  value  should  approach  unity;  but  especially  with  the  advent  of 


supercomputers  with  novel  architecture,  or  with  I/O  transfer  to  external  memory  devices.  It  could  well 
be  2.0  or  even  larger. 
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The  number  of  arithmetic  operations  per  grid  point  per  Iteration,  Wgj,  Is  a  strong  function  of  the 
numerical  method;  that  Is,  of  the  flow  equations,  the  boundary  conditions,  the  solution  algorithm,  and 
whether  the  grid,  metrics  are  computed  at  each  time  step  or  stored  In  memory.  The  quantity  Ng  represents 
the  number  of  grid  points  for  a  finite-difference  method,  the  number  of  elements  for  a  finite-element 
method,  or  the  number  of  panels  for  a  panel  method.  Consequently,  Wg]  Ng  represents  the  number  of  arith¬ 
metic  operations  that  must  be  performed  at  each  Iteration  or  time  step,  although  In  some  Instances  with 
panel  methods,  Ng  log  Ng  Is  a  more  accurate  representation  than  N*. 

t  The  total  number  of  Iterations,  N j ,  to  converge  to  the  desired  level  of  accuracy  (or  the  total  num¬ 
ber  of  time  steps  In  an  unsteady  problem).  Is  generally  the  most  difficult  quantity  to  estimate.  And 
unfortunately,  N[  for  a  given  algorithm  Is  often  both  problem-dependent  and  grid-dependent.  This  Is 
particularly  true  In  the  early  stages  of  maturtty  of  a  code,  such  as  many  of  those  presented  In  Sessions 
III  -  V.  The  value  of  N[  Is  generally  significantly  greater  for  explicit  methods  than  for  Implicit 
ones;  but  on  the  other  hand,  Wgj  for  Implicit  methods  tends  to  be  greater. 

Finally,  the  sustained  computing  speed,  FLOPS,  Is  a  function  of  the  computer  clock  speed  and  archi¬ 
tecture,  the  data  structure  and  techniques  of  the  code,  and  the  memory  requirements  (In-core  or  external 
memory).  The  aerodynamlclst  has  relatively  little  control  over  this  quantity;  rather,  It  Is  largely  dic¬ 
tated  by  management,  l.e.,  by  the  size  and  cost  of  the  computer  system. 

It  Is  clear  from  Eqn.  1  that  many  different  factors  determine  the  CPU  time,  and  hence,  the  computa¬ 
tional  cost,  of  an  aerodynamic  calculation.  Unfortunately,  none  of  the  papers  provided  all  of  the  Ingre¬ 
dients  specified  In  Eqn.  1;  however,  enough  Information  was  provided  e*ther  at  the  Symposium  or  elsewhere, 
or  could  be  Inferred,  to  establish  some  trends  and  order-of-magnltude  estimates.  This  Information  Is 
summarized  In  Table  1,  for  a  hypothetical  wing-body  combination  of  moderate  geometrical  complexity.  Here 
the  Mach  number  Is  Implicitly  assumed  to  be  subsonic  In  the  linear  case,  supersonic  for  the  Parabolized 
Navler-Stokes  estimate,  and  transonic  otherwise;  and  the  reference  computer  speed  Is  that  of  the  Cray  1-S, 
l.e.  FLOPS  =  80  »  106. 

It  should  be  emphasized  that  these  estimates  are  very  approximate,  hypothetical,  and  somewhat 
abrltrary;  therefore,  they  could  easily  be  off  by  a  factor  of  2  or  more.  Also,  they  reflect  the  arguable 
premise  that  Increasing  numbers  of  grid  points  should  be  accompany  the  Increasing  sophistication  In  the 
flow  modeling.  In  order  to  capture  the  more  complex  flow  physics  that  motivate  the  more  complex 
approaches.  Nevertheless,  Table  1  gives  a  general,  qualitative  picture  of  the  differences  In  the  overall 

Table  l.  APPROXIMATE  COMPUTATIONAL  REQUIREMENTS 
FOR  A  COMPLEX  WING-BOOT  COMBINATION 


Method 

«GI 

V 

millions 

NI 

CPU, 

minutes 

Total  memory, 
words  xlO6 

Linear  (panel) 

(1000  panels) 

2  -  20 

0.5  -  1.0 

Transonic  Small 
Disturbance 

100 

0. 1-0.2 

100  -  300 

5  -  15 

0.5  -  1.0 

Full  Potential 

500 

0.1-0. 2 

200  -  500 

10  -  30 

1  -  4 

Euler 

1000- 

3000 

0.2-0. 5 

500  -  5000 

50  -  500 

2  -  10 

Parabolized 

Navler-Stokes 

Ng  x  N,  -  2-5xl06 

10  -  60 

0.5  -  1.0 

Navler-Stokes 

1500- 

4500 

0.5-2 

1500  -  10000 

1000  -  5000 

15  -  60 

computational  requirements  of  many  of  the  approaches  that  were  described  at  the  Symposium.  It  Indicates, 
for  example,  that  transonic  small-disturbance,  full-potential,  and  PNS  codes  can  be  more  or  less  competi¬ 
tive  with  panel  codes,  with  respect  to  computer  resources.  The  table  also  gives  some  Idea  of  the  price 
one  would  have  to  pay  today  to  obtain  full  viscous  simulations  of  transonic  flow.  These  observations 
deserve  further  conment,  and  this  leads  logically  to  the  question  of  what  capabilities  the  different 
methods  actually  offer  for  aerodynamic  design  and  analysis. 

3.1.2  Capabilities  and  Limitations.  Linear  panel  codes  have  been  the  workhorses  of  the  cosputa- 
tlonal  aerodynamics  cormunlty  for  many  years,  and  they  are  generally  considered  to  be  the  most  flexible 
and  least  Inexpensive  approach  for  complex  configurations.  They  are  Inappropriate,  however,  for  transonic 
flows;  and  as  noted  above,  nonlinear  potential  codes  have  been  developed  to  a  competitive  position 
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costwlse.  Additional  difficulties  can  arise  In  supersonic  flows,  cf.  Paper  No.  5;  therefore,  the 
estimates  of  computer  resources  In  Table  1  suggest  that  space-marching  Euler  or  PNS  codes  might  well  be 
attractive  alternatives  for  supersonic  cases.  However,  these  nonlinear  codes  have  had  less  time  to 
benefit  from  user  experience,  and  they  require  more  effort  to  generate  the  grids.  Consequently,  they 
would  probably  be  more  difficult  to  use  today  In  an  engineering  environment. 

Today,  almost  any  company,  organization,  or  Institution  that  has  significant  connections  with  tran¬ 
sonic  aircraft  or  missiles  routinely  uses  nonlinear  3-0  potential  codes,  and  the  helicopter  comnunlty  is 
moving  rapidly  In  that  direction.  The  recent  trend  to  move  on  to  Euler  codes  seem  to  be  motivated  by 
three  main  factors:  first,  to  treat  stronger  shock  waves  than  Is  permitted  by  the  potential  formulation; 
second,  the  ability  to  capture  vortex  sheets  and  other  aspects  of  rotational  flows;  and  third,  to  build  a 
bridge  to  Navler-Stokes  codes.  The  results  of  AGARO  Working  Group  WG-07  [ Ref.  7]  are  especially  note¬ 
worthy  with  regard  to  the  first  consideration,  and  Papers  14,  15,  and  16,  as  well  as  Ref.  7,  provide 
representative  demonstrations  of  the  second  factor. 

Two  points  can  be  made  regarding  the  nonlinear  Invlseld  codes,  especially  the  Euler  ones.  With  the 
neglect  of  viscous  effects,  the  calculated  lift  Is  too  high,  the  drag  Is  too  tow,  and  the  shock  wave  Is  In 
the  wrong  position;  consequently,  the  pitching  moment  Is  generally  Incorrect,  too.  In  practice,  these 
effects  may  be  masked  or  counterbalanced  by  the  effects  of  coarse  grids,  numerical  dissipation,  and  other 
numerical  errors  (as  Illuminated  by  Paper  14,  for  example).  The  net  effect  Is  to  Introduce  an  element  of 
uncertainty  In  the  results,  which  can  seldom  be  evaluated  short  of  performing  extensive,  and  expensive, 
grid-refinement  studies.  Similar  criticisms  can  be  leveled  at  the  Navler-Stokes  calculations  as  well,  as 
discussed  below. 

Secondly,  the  rise  In  applications  to  complex  3-D  configurations  has  dramatically  Increased  the 
Importance  of  both  grid  generation  and  post-processing  data  analysis,  to  such  an  extent  that  now  they  ore 
often  the  pacing  Items  with  respect  to  manpower  and  elapsed  time.  As  Dr.  Wolfgang  Schmidt  noted  In  the 
closing  discussion  period.  It  may  take  months  to  set  up  the  complex  mesh  for  a  complete  aircraft,  followed 
by  a  day  or  two  of  “clock”  time  to  complete  an  hour  or  so  of  CPU  time,  and  It  may  then  take  months  to 
“reduce"  the  data  fully.  (The  comparison  with  large  wind-tunnel  projects  goes  without  saying.)  Further¬ 
more,  this  added  complexity  means  that  It  Is  becoming  Increasingly  difficult  to  hand-off  these  powerful 
nonlinear  codes  to  Inexperienced  users. 

Coupled  vlscous-lnvtscld  Interaction  methods  are  generally  considered  to  be  promising  for  design 
applications  because  of  their  computational  efficiency  relative  to  Navler-Stokes  approaches.  Two- 
dimensional  formulations  linking  Invlseld  and  viscous  algorithms  formed  the  basis  of  a  major  AGARD  Sympo¬ 
sium  In  1980  |Ref.  8|.  These  methods  have  matured  and  come  Into  general  use  since  then.  Including 
successful  applications  to  flows  with  small  amounts  of  separation.  Several  successful  extensions  to  3-D 
have  appeared  In  the  recent  literature,  although  these  are  mostly  weak-coupllng  methods  applied  to  unsepa¬ 
rated  flows.  However,  the  confidence  level  that  exists  for  the  2-D  methods  seemed  to  be  lacking  at  this 
Symposium.  As  noted  above,  this  may  reflect  a  period  of  renewed  effort  to  develop  new  strong-coupling 
methods  that  are  robust  and  efficient,  and  as  noted  by  Fermln  In  Paper  26,  which  give  accurate  predictions 
of  drag.  The  ensemble  of  papers  In  Session  IV  suggests  that  much  work  remains  to  be  done,  even  for  high- 
aspect-ratio  wings,  when  strong  coupling  Is  required  between  the  viscous  and  Invlseld  regions.  The  new 
RAE  experiments  discussed  In  Paper  26  should  provide  good  targets  and  challenges  for  Improvements  In  this 
area. 

Within  the  scope  of  tne  Reynolds-averaged  formulation  of  the  Navler-Stokes  equations,  the  limitations 
of  the  preceding  methods  are  theoretically  eliminated.  The  computational  requirements  shown  in  Table  1 
Indicate  that  this  approach  Is  not  yet  practical  for  routine  aeronautical  analysis  and  design,  although 
this  conclusion  would  completely  change  within  a  decade  If  the  past  trends  In  computer  technology  and 
algorithm  development  are  extrapolated  Into  the  future  ( Ref.  9).  However,  the  primary  Issue  here  Is  the 
validity  of  the  results  that  are  being  computed  today.  In  most  cases,  2-D  Navler-Stokes  calculations  are 
approximately  as  good  as  the  average  wind-tunnel  results  (Ref.  101.  However,  the  Issue  of  3-0  code  vali¬ 
dation  Is  mired  In  controversy  and  uncertainty  regarding  turbulence  modeling,  spatial  resolution,  the 
effects  of  numerical  dissipation  and  other  numerical  errors,  and  the  completeness  and  reliability  of  the 
relevant  experiments. 

For  most  relatively  benign  flows,  the  Navler-Stokes  calculations  agree  reasonably  well  with  experi¬ 
ments  (as  do  the  results  of  the  simpler  codes  at  a  fraction  of  the  cost),  provided  “reasonable"  grids  are 
used  and  some  attention  Is  paid  to  the  numerical  dissipation  parameters.  However,  In  those  cases  where 
the  Navler-Stokes  approach  Is  fully  Justified  over  simpler  methods  because  of  strong  vlscous-lnvlscld 
Interaction  and/or  massive  flow  separation,  virtually  all  of  results  to  date  have  showed  significant  dis¬ 
crepancies  when  detailed  comparisons  were  made  with  experiments.  In  the  opinion  of  most  turbulence 
experts,  the  simple  turbulence  models  used  in  all  the  large  Navler-Stokes  codes  described  at  the  Symposium 
are  Inadequate  for  such  applications,  and  In  the  opinion  of  the  leading  numerical  analysists,  the  grids 
would  have  to  be  refined  considerably  before  the  solutions  could  be  expected  to  be  grid-independent. 

Therefore,  one  may  conclude  that  today's  3-D  Navler-Stokes  solutions  are  useful  but  largely  qualita¬ 
tive.  However,  the  extent  of  the  Information  burled  In  the  results  Is  enormous,  and  they  probably  contain 
valuable  physical  Insights  that  may  not  be  evident  In  experimental  data.  The  Issues  of  turbulence  model¬ 
ing  and  grid  resolution  are.  In  fact,  being  addressed  vigorously,  and  the  papers  In  Session  V  give  a  pre¬ 
view  of  some  things  expect  In  the  future. 
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3.2  Additional  Problem  and  Issues 

In  Addition  to  the  Issues  of  costs  and  capabilities  of  the  various  CFD  approaches,  the  topics  listed 
below  emerged  either  during  the  formal  sessions  or  In  the  final  discussion  period. 

3.2.1.  Adaptive  Grids.  None  of  the  regular  papers  discussed  solution-adaptive  grids.  However,  this 
concept  offers  the  potential  of  reducing  the  computational  requirements  by  placing  grid  points  where  they 
are  most  needed  and  making  more  efficient  use  of  the  computer  resources  available.  Several  participants 
Indicated  that  this  topic  Is  under  active  Investigation. 

3.2.2.  Turbulence  Model  vs.  Numerical  Errors.  When  numerical  solutions  are  suspected  to  be  In 
error,  one  may  question  either  the  turbulence  model  or  the  numerical  model.  This  Issue  has  several  conso¬ 
nants.  The  first  Is  the  relative  magnitude  of  the  numerical  viscosity  that  Is  Inherent  In  most  CFD 
methods  compared  to  the  laminar  or  turbulent  viscosity;  one  may  ask  which  of  these  Is  the  dominant  factor 
In  different  parts  of  the  flow  field.  The  numerical  viscosity  normally  varies  In  proportion  to  the  grid 
spacing,  and  thus  Is  a  property  of  the  numerics,  whereas  the  physical  viscosity  Is  a  property  of  the  fluid 
and/or  flow  gradients.  In  the  Euler  calculations  of  Murman  and  Rlzzl,  for  example,  the  numerical  viscos¬ 
ity  Is  the  mechanism  which  permits  the  formation  of  shear  layers  In  a  numerical  simulation  of  an  Invlscld 
flow.  However,  they  claimed  that  any  amount  of  numerical  viscosity,  however  small,  was  sufficient  to 
allow  a  realistic  solution  to  develop,  and  that  their  results  were  Insensitive  to  Its  value.  In  some  of 
the  Navler-Stokes  calculations,  (Papers  27  and  30,  for  example),  the  numerical  viscosity  was  rather  large 
In  some  parts  of  the  flow  field,  especially  near  shock  waves,  but  was  probably  small  compared  to  the  tur¬ 
bulent  viscosity  In  the  interior  of  the  boundary  layer.  Unfortunately,  It  has  not  yet  been  possible  to 
refine  these  grids  enough  to  determine  how  sensitive  these  3-D  solutions  are  to  numerical  viscosity,  or 
other  numerical  parameters. 

Another  Issue,  raised  by  Or.  Tuncer  Cebecl,  Is  the  grid  resolution  In  the  direction  normal  to  the 
wall  that  Is  required  for  accurate  prediction  of  the  surface  shear  stress.  Boundary  layer  methods  tend  to 
use  many  more  points  across  the  viscous  layer  than  Navler-Stokes  methods,  although  the  trends  are  to  use 
comparable  spacing  Immediately  adjacent  to  the  wall.  Again,  more  grid-refinement  studies  are  needed. 

With  the  exception  of  Papers  19,  25,  and  26,  all  of  the  turbulent  calculations  were  done  using  eddy- 
viscosity  models  for  simplicity.  As  noted  In  Section  3.1.2,  this  would  be  a  source  of  serious  error  for 
many  complex  flows.  In  the  opinion  of  most  turbulence  experts.  However,  this  Issue  was  not  a  major  theme 
of  the  Symposium,  nor  an  overriding  concern  of  most  of  the  participants. 

3.2.3.  Drag  Calculations.  This  topic  was  given  only  secondary  treatment  In  almost  all  of  the  formal 
papers.  However,  In  the  discussion  period  the  accuracy  of  drag  predictions  emerged  as  being  very  Impor¬ 
tant  to  many  participants,  and  It  seems  certain  to  receive  more  attention  In  the  future.  Mr.  J.  w.  Slooff 
stated,  with  reference  to  the  large  number  of  colored  graphs  of  results  that  were  presented,  that  an 
Improvement  In  drag  prediction  of  10  counts  would  be  worth  more  than  a  thousand  pictures. 

3.2.4.  Hlqhly-Vortlcal  Flows.  A  major  motivation  for  many  of  the  advanced  CFD  methods  under  devel¬ 
opment  today  Is  the  desire  to  compute  the  entire  flow  field  of  flight  vehicles  at  high  angles  of  attack, 
and  these  flows  contain  complex  vortical  structures.  Such  flow  features  present  major  challenges  In  grid 
generation  and  computation,  as  many  existing  techniques  are  not  suitable  for  capturing  and  preserving  the 
strong  flow  gradients  that  are  Involved.  Solution-adaptive  grid  techniques  may  be  advantageous  here. 

Also,  new  techniques  of  capturing  or  fitting  concentrated  vortices  may  have  to  be  developed,  by  analogy 
with  the  successful  efforts  In  the  past  to  develop  both  shock-fitting  and  shock-capturing  methods. 

This  class  of  flows  also  provides  major  challenges  In  analyzing  and  understanding  the  computed 
results,  as  the  topologies  of  vortical  structures  associated  with  3-D  separation  are  very  complex. 

Paper  31  Included  some  exciting  examples  of  this  complexity,  which  requires  extensive  computer-graphics 
capability  to  even  begin  to  analyze  the  flow  structure. 

3.2.5.  Future  Topics.  The  fixed-wing  aircraft  Industry  has  been  the  primary  driving  force  for  the 
computational  aerodynamics  techniques  that  were  presented  at  this  Symposium.  As  CFD  methods  for  this 
class  of  vehicles  mature,  one  may  expect  to  see  expanded  applications  to  rotorcraft  and  v/STOl  aircraft, 
turbomachinery,  and  hypersonic  vehicles.  Participants  also  suggested  that  major  future  efforts  may  Include 
transition  prediction;  unsteady  aerodynamic  flows;  advanced  turbulence  modeling,  Including  numerical  simu¬ 
lations  of  turbulence;  and  optimization  of  body  shapes  to  satisfy  prescribed  aerodynamic  properties,  such 
as  pressure  distributions. 

Finally,  In  Paper  27,  Dr.  J.S.  Shang  discussed  at  length  "Interdisciplinary  computational  fluid 
dynamics,”  In  which  fluid-flow  equations  (Navler-Stokes  equations,  In  his  perspective)  are  coupled  with 
the  governing  equations  of  solid  mechanics,  chemistry,  combustion,  electromagnetics,  optics,  etc. 

4.  CONCLUSIONS  AND  RECOMMENDATIONS 

Computational  Fluid  Oynimlcs  has  established  a  firm  role  In  aerodynamic  design,  and  this  Symposium 
revealed  the  wide  spectrum  of  CFO  tools  that  are  now  available  to  the  aeronautics  community.  Furthermore, 
the  application  of  these  tools  to  problems  of  ever-increasing  complexity  Is  accelerating  rapidly,  and 
computer  hardware  Is  advancing  even  faster  than  Improvements  In  numerical  algorithms.  As  a  result,  the 
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users  of  modern  CFD  technology  can  generate  vast  amounts  of  Information,  and  their  managers  are  not  always 
sure  how  valid  or  useful  this  mountain  of  data  really  Is. 

Three  Important  needs  arise  from  this  situation.  The  first  Is  to  raise  the  confidence  level  of  CFO 
results,  as  opposed  to  overselling  their  Importance  or  covering  up  their  limitations.  That  Is,  we 
urgently  need  better  validation  of  the  codes  and  convincing  demonstrations  of  their  capabilities  and 
limitations.  This  requires  more  careful  grid-ref Inement  studies  and  numerical -error  checks,  leading  to 
the  establishment  of  reliable  error  bands  on  CFO  results.  Uncertainty  analysis  has  become  an  Important 
part  of  quality  experimental  research  and  testing,  and  there  Is  no  reason  not  to  establish  similar  quality 
control  for  CFO.  It  Is  also  essential  to  establish  more  standard  test  cases  and  to  compare  results  for 
them.  A  limited  number  of  new  experimental  configurations  should  be. established  with  the  specific  needs 
of  CFO  code  validation  In  mind,  and  with  more  redundancy,  higher  accuracy,  and  better  flow  quality  than 
has  been  the  norm  heretofore.  Numerous  AGARD  working  groups  and  other  parties  have  performed  valuable 
services  of  this  type  In  recent  years,  by  defining  test  cases,  coordinating  the  experimental  efforts  at 
several  different  laboratories,  and  assessing  the  results.  This  kind  of  effort  should  be  continued  with 
renewed  vigor. 

Second,  better  data  and  Information  management  tools  are  becoming  as  Important  as  better  computing 
hardware  for  post-processing  the  CFO  results  and  extracting  the  valuable  Information  that  they  contain. 
Although  not  all  of  the  Symposium  participants  shared  the  writer's  enthusiasm  for  color  computer  graphics, 
they  all  agreed  that  this  aspect  of  applied  CFO  Is  essential  and  Is  growing  In  Importance.  This  Is 
another  area  In  which  AGARO  should  try  to  facilitate  better  multi-national  cooperation  and  Information 
exchange. 

Third,  the  Issue  of  turbulence  modeling  continues  to  be  a  cloud  over  CFD  results  for  viscous  prob¬ 
lems.  Although  this  Issue  was  not  a  major  theme  of  this  Symposium,  It  will  probably  hinder  both  progress 
In  and  acceptance  of  applied  computational  aerodynamics  for  many  years.  On  the  one  hand,  CFD  researchers 
must  forge  ahead  In  the  development  of  new  methods,  without  waiting  for  the  fundamentals  of  turbulence  to 
be  understood  fully.  However,  they  need  to  structure  their  codes  as  much  as  possible  to  accept  new  turbu¬ 
lence  models,  and  they  must  be  willing  to  Incorporate  them  as  better  ones  become  available. 

In  conclusion,  the  Symposium  provided  a  timely  and  valuable  forum  for  exchanging  Information  about 
recent  developments  In  applied  computational  aerodynamics,  and  It  clearly  fulfilled  the  stated  goal  of 
assessing  the  status  of  CFO  In  contemporary  aeronautical  design  and  analysis.  Together  with  the  problems 
and  Issues  were  that  were  Identified,  many  successes  of  CFD  In  recent  years  were  highlighted,  and  the 
trends  of  current  developments  point  to  a  promising  future. 

In  the  near  term,  refinements  In  the  vlscous-lnvlscld  methods  should  overcome  many  of  the  shortcom¬ 
ings  of  present  Invtscld  methods  without  significant  additional  costs,  and  continuing  Improvements  In 
computer  hardware  and  software  and  In  numerical  algorithms  will  create  more  demand  for  the  Navler-Stokes 
codes.  Thus  the  present  level  of  maturity  of  Invlscld  computational  aerodynamics  can  be  expected  to 
extend  to  viscous  flows  within  a  few  years  for  fixed-wing  aircraft.  At  the  same  time,  CFD  should  find 
rapidly-increasing  applications  In  other  fields  of  aeronautics,  such  as  rotorcraft  and  turbomachinery.  In 
the  longer  term.  It  seems  clear  that  CFO  will  be  combined  with  a  wide  range  of  other  disciplines  to  expand 
greatly  Its  range  of  applications  and  usefulness. 
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